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Using  numerical  modeling,  it  is  shown  that  chains  of  dielectric  spheres  and  cylinders  act  as  polarizers.  The  mecha¬ 
nism  is  based  on  gradual  filtering  of  periodically  focused  modes  with  a  certain  polarization  propagating  with  min¬ 
imal  losses  due  to  Brewster  angles  conditions,  whereas  orthogonally  polarized  modes  are  strongly  attenuated.  It  is 
shown  that  chains  of  cylinders  filter  linearly  polarized  beams,  whereas  chains  of  spheres  filter  radially  polarized 
beams.  In  the  geometrical  optics  limit,  we  show  that  in  a  range  of  sphere  refractive  indices  1.68-1.80  a  degree  of 
radial  polarization  in  excess  of  0.9  can  be  obtained  in  10-sphere-long  chains.  ©  2013  Optical  Society  of  America 
OCIS  codes:  (120.3150)  Integrating  spheres;  (130.5440)  Polarization-selective  devices;  (350.3950)  Micro-optics. 
http://dx.doi.Org/10.1364/OL.38.004208 


Recently,  it  was  demonstrated  that  linear  arrays  of  di¬ 
electric  spheres  can  operate  as  lossless  waveguides  for 
beams  with  certain  spatial  and  polarization  properties 
[1].  This  can  be  achieved  due  to  the  fact  that  the  Brewster 
angles  conditions  are  periodically  reproduced  for  modes 
with  a  period  matching  the  size  of  two  sphere  diameters 
{2D).  In  the  geometrical  optics  limit,  it  was  shown  that 
such  periodically  focused  modes  (PFMs)  are  radially 
polarized. 

It  should  be  noted  that  the  PFM  concept  has  been 
stimulated  by  previous  studies  of  the  chains  of  mesoscale 
(D  <  lOi,  where  X  is  the  wavelength  of  light)  spheres  [2] 
and  cylinders  [3]  performed  in  the  context  of  coupling 
between  their  whispering  gallery  modes,  beam  focusing 
produced  by  individual  spheres  and  microdisks  [4-7] 
(termed  photonic  nanojet  effect),  and  coupling  between 
nanojets  in  chains  of  spheres  [8].  The  latter  mechanism 
leads  to  formation  of  periodic  nanojet-induced  modes 
(NIMs)  observed  in  chains  of  polystyrene  microspheres 
[9-13]  and  cylinders  [14].  In  principle,  NIMs  in  mesoscale 
spheres  can  be  considered  as  a  wave-optics  analog  of 
PFMs  introduced  [1]  in  a  geometrical  optics  limit.  Some 
common  NIMs  and  PFMs  properties  leading  to  a  gradual 
reduction  of  the  focused  beam  width  [1,9]  along  the  chain 
and  a  gradual  reduction  of  the  attenuation  below 
0.1  dB/sphere  [1,1^]  have  been  reported. 

In  this  work  we  show  that  chains  of  cylinders  and 
spheres  can  also  be  used  as  polarizers  with  focusing 
capability.  The  operation  of  such  structures  is  based 
on  the  fact  that  their  photonic  eigenstates  are  repre¬ 
sented  by  PFMs  with  orthogonal  polarizations,  which 
have  dramatically  different  propagation  losses.  We  first 
illustrate  these  effects  by  finite  element  method  (FEM) 
simulations  for  coupled  cylinders  where  the  eigenstates 
are  linearly  polarized.  After  that  we  consider  these  ef¬ 
fects  in  coupled  microspheres  where  the  polarization  ei¬ 
genstates  are  represented  by  radially  and  azimuthally 
polarized  beams.  Using  ray  tracing  for  sufficiently  large 


spheres  {D  »  lOi)  we  show  that  PFMs  with  small  losses 
can  be  obtained  in  a  broad  range  of  refractive  indices 
from  ^/2  to  2.  We  show  that  high  degrees  of  polarization 
(DOP),  --0.9,  can  be  obtained  in  10-sphere  long  chains 
using  different  ray  sources.  The  polarization  filtering 
and  focusing  properties  display  an  interesting  interplay 
in  chains  of  spheres  so  that  higher  degrees  of  radial 
polarization  favor  sharper  focus  [1^,1^]  with  applications 
in  laser  scalpels  [17,1^],  optical  tweezers  [1^],  and  vari¬ 
ous  focusing  devices  [20,21]. 

The  results  of  two-dimensional  (2D)  FEM  modeling  of 
polarization-dependent  transmission  through  a  series 
of  dielectric  cylinders  with  refractive  index  n  =  Vs 
performed  by  COMSOL  Multiphysics  are  illustrated  in 
Fig.  1.  Launching  plane  waves  polarized  along  the  y  axis 
leads  to  higher  transmission  compared  to  the  case  of 
incident  waves  polarized  along  the  ^  axis,  as  can  be  seen 
by  comparing  Figs.  1(a)  and  1(b),  respectively.  In  both 
cases  the  typical  PFM  pattern  [1]  is  observed.  Higher 
transmission  for  incident  polarization  along  the  y  axis 
can  be  explained  by  the  role  of  the  Brewster  angle  for 
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Fig.  1.  2D  electric  field  intensity  maps  calculated  for  touching 
cylinders  with  D  =  3  |im  and  n  =  a/3  at  A  =  0.53  |im  for  the 
incident  plane  waves  with  electric  field  oscillating  along 
the  (a)  y  axis  and  (b)  axis,  (c)  The  intensity  profiles  of 
the  focused  beams  calculated  0.62  pm  away  from  the  tip  of  the 
end  cylinder  shovring  stronger  peak  transmission  in  case 
(a)  compared  to  (b). 


0146-9592/13/204208-04$15.00/0 


©  2013  Optical  Society  of  America 


October  15,  2013  /  Vol.  38,  No.  20  /  OPTICS  LETTERS  4209 


this  polarization.  Since  only  a  small  fraction  of  the  wave- 
front  is  incident  on  the  cylinder  surface  at  the  Brewster 
angle,  this  model  cannot  be  used  for  quantitative  analysis 
of  such  mesoscale  structures;  however  it  provides  a 
physical  explanation  of  this  observation. 

Similar  ideas  can  be  applied  in  a  more  elaborate  way 
to  structures  formed  by  macroscopic  spheres  {D  »  lOi) 
using  geometrical  optics  approximation.  Previously, 
we  showed  that  the  exact  conditions  for  Brewster  inci¬ 
dence  [1]  can  be  satisfied  for  spheres  with  refractive 
index  n  =  V3  if  collimated  TM  polarized  incident  rays 
are  off-axially  shifted  by  r/D  =  V3/4,  as  illustrated  in 
Fig.  2(a).  It  is  interesting,  however,  to  consider  different 
spheres’  indices  and  study  if  PFMs  can  still  exist  under 
these  conditions.  The  idea  of  our  analysis  is  that  in  order 
to  remain  2D  periodic,  such  modes  should  have  different 
off-axial  shifts.  Their  external  {6i)  and  internal  an¬ 
gles  of  incidence  and  refraction,  respectively,  would 
not  be  exactly  equal  to  those  following  Brewster’s  law, 
but  they  can  be  sufficiently  close  to  the  optimal  values. 
Most  importantly,  these  angles  would  be  periodically  re¬ 
produced  as  light  propagates  through  the  chains  which 
would  help  to  keep  the  losses  low. 

The  collimated  incident  rays  are  reproduced  with 
2D  period  when  the  refracted  ray  passes  through  the 
point  where  the  spheres  touch  leading  to  the  following 
equations:  6i  =  2  cos“^(?^/2?^^),  r/D  =  0.5  sin(^^),  where 
=  1  is  the  refractive  index  of  the  background  medium. 
Based  on  these  equations,  it  can  be  deduced  that 
quasi-PFMs  exist  in  a  broad  range  (\/2  <  n  <  2)  of  indi¬ 
ces.  The  angle  of  incidence  and  axial  offset  of  these 
modes  decrease  with  the  increase  in  the  refractive  index. 

The  attenuation  properties  of  chains  of  spheres  were 
studied  by  considering  rays  forming  configurations  with 
2D  period.  For  each  index  it  was  achieved  by  adjusting  6i 
and  r  according  to  the  equations  presented  above.  Propa¬ 
gation  losses  for  the  incident  rays  illustrated  in  Fig.  2(a) 
can  be  estimated  by  taking  into  account  Fresnel  reflec- 


Fig.  2.  (a)  Ray  tracing  in  an  array  of  identical  touching  spheres 

under  the  condition  of  periodic  propagation,  (b)  Transmittance 
and  (c)  loss  in  chains  of  =  1,  5,  10  and  20  spheres  as  a  func¬ 
tion  of  n  for  TE  and  TM  polarizations  of  the  incident  ray.  In  a 
geometrical  optics  limit  (D  »  lOA)  the  results  are  not  depen¬ 
dent  on  D  and  X. 


tion  coefficients  [^]  at  each  spherical  interface  from 
the  following  equations  for  TE  and  TM  polarizations: 


sin2(36>^/2)  ’ 


(1) 


_  tan^(6>^/2) 

“  tan2(36>^/2)  ’ 


(2) 


The  transmittance  and  corresponding  reflection  losses 
for  rays  with  TE  and  TM  polarizations  as  a  function  of  n 
for  chains  ofN  =  1,  5, 10,  and  20  spheres  are  presented  in 
Figs.  2(b)  and  2(c),  respectively.  As  expected,  regardless 
of  the  length  of  the  chain,  full  transmission  and  zero  re¬ 
flection  loss  occurs  at  n  =  VS  for  TM  polarization  which 
satisfies  the  Brewster  angles  conditions  [1].  It  is  interest¬ 
ing,  however,  that  the  high  transmission  properties  are 
preserved  in  a  broad  range  of  indices  around  n  = 
for  TM  polarization.  It  is  illustrated  in  Fig.  2(c)  that  for 
10-sphere  long  chains  with  1.68  <  n  <  1.80,  TM  polarized 
PFMs  have  total  propagation  losses  smaller  than  1  dB, 
i.e.,  less  than  0.1  dB/sphere.  Conversely,  for  TE  polariza¬ 
tion,  the  total  PFM  losses  exceed  20  dB  in  the  same  range 
of  indices.  This  means  that  for  randomly  polarized  colli¬ 
mated  incident  rays,  the  output  mainly  has  a  TM  polari¬ 
zation  component.  Due  to  axial  symmetry  of  the  problem, 
the  global  state  of  polarization  of  the  transmitted  beam 
should  be  increasingly  radial. 

To  study  the  polarizing  capability  of  such  chains  we 
considered  three  types  of  ray  sources: 

(i)  A  source  of  collimated  rays  with  uniform  density 
across  the  structure  and  random  polarization,  as  shown 
in  Fig.  3(a).  Such  rays  impinge  on  the  spherical  surface  at 
different  angles  depending  on  r.  Only  the  TM  polarized 
rays  with  parameters  6i  and  r  close  to  those  presented 
in  the  above  equations  can  effectively  contribute  to 
the  transmission  through  sufficiently  long  chains.  Such 
source  effectively  represents  a  geometrical  optics  model 
for  illumination  with  a  collimated  laser  beam. 

(ii)  A  sphere  with  diameter  D  that  emits  rays  through¬ 
out  its  volume  in  all  directions  with  random  polariza¬ 
tion  as  shown  in  Fig.  3(b).  This  light  source  (S)  can 
be  considered  as  a  model  for  dye-doped  fluorescent 


Fig.  3.  Three  different  types  of  coupling  to  the  PFMs. 
(a)  Collimated  rays,  (b)  spherical  emitter,  and  (c)  a  multimode 
fiber  inserted  inside  a  hollow  waveguide. 
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microspheres  which  have  been  used  in  experimental 
studies  [1,9-1^]  of  such  structures  below  the  lasing 
threshold  for  whispering  gallery  modes.  Above  this 
threshold  the  geometrical  optics  approximation  typically 
cannot  be  used  and  full-wave  EM  modeling  is  re¬ 
quired  [3,^]. 

(iii)  A  source  of  diverging  rays  mimicking  the  direction¬ 
ality  of  multimode  fibers  with  typical  divergence  angle 
around  12°.  Such  directionality  can  be  realized,  as  an  ex¬ 
ample,  in  the  mid-infrared  range  using  fibers  with 
'-ISO  pm  core  diameters,  as  shown  in  Fig.  3(c).  This  type 
of  illumination  can  take  place  in  focusing  devices  oper¬ 
ating  with  multimodal  beams  [17,1^].  This  light  source 
represents  intermediate  directional  properties  between 
cases  (i)  and  (ii). 

The  DOP  P(r)  is  defined  as  the  ratio  of  the  (averaged) 
intensity  of  the  polarized  portion  of  the  beam  to  its  total 
(averaged)  intensity,  both  taken  at  the  same  point  [M]. 
This  rigorous  and  unambiguous  definition  of  the  DOP 
is,  however,  sometimes  difficult  to  use  in  practical  cases. 
In  such  situations  various  ad  hoc  definitions  of  the  DOP 
are  frequently  used  in  the  form  [M]: 


|4(r)  -4(r)| 


4(r)  +  lyir)  ’ 


(3) 


where  and  ly  are  the  averaged  intensities  in  two  mu¬ 
tually  orthogonal  directions,  their  choice  being  suggested 
by  the  geometry  of  the  problem.  Unlike  the  P(r),  the 
quantity  0(r)  depends  on  the  choice  of  the  x,  y  axes. 
It  has  been  theoretically  demonstrated,  however,  that 
the  ad  hoc  definition  in  the  form  of  Eq.  (^  will  correctly 
represent  the  degree  of  polarization  P(r),  provided  that 

and  ly  are  taken  to  be  true  eigenvalues  of  the  polari¬ 
zation  matrix  representing  the  symmetry  of  a  given  prob¬ 
lem  [M].  The  axial  symmetry  of  the  chains  of  spheres 
implies  that  their  polarization  eigenvalues  should  be  rep¬ 
resented  by  the  intensities  of  radially  (/^)  and  azimuthally 
(/^)  polarized  beams. 

Typical  intensity  distribution  calculated  for  spheres 
with  n  =  Vs  using  a  flat  detector  in  contact  with  the 
end  sphere  is  illustrated  in  Fig.  4(a)  in  the  case  of  colli¬ 
mated  incident  rays  for  20-sphere  long  chain.  This  case  is 
representative  of  the  calculations  performed  for  chains 
formed  by  even  number  of  spheres  with  n  =  1.65-1.85. 
For  this  range  of  indices  the  output  beams  consist  of  a 
ring  with  radius  r  determined  by  the  PFM  axial  offset, 
a  more  compact  central  beam  formed  by  paraxial  rays, 
and  a  weak  background  illumination.  The  dominant  con¬ 
tribution  to  the  output  optical  power  is  given  by  the  ring 
with  radius  r. 

To  estimate  the  DOP,  we  placed  a  horizontally  oriented 
linear  polarizer  (along  x')  after  the  end-sphere,  as  shown 
in  Fig.  4(b).  The  azimuthal  intensity  modulation  along  the 
ring  in  Fig.  4(c)  carries  information  about  DOP  of  the 
transmitted  beam.  It  was  studied  by  using  a  local  detector 
(extended  along  x')  with  dimensions  D/lOxD/lOO 
placed  after  the  polarizer,  as  schematically  indicated 
by  a  stripe  labeled  in  Figs.  4(b)  and  4(c).  Such  a  de¬ 
tector  provides  a  local  calculation  of  the  intensity  of 
TM  components  of  transmitted  PFMs,  which  in  turn 
can  be  taken  as  a  measure  of  From  the  axial  symmetry 


Fig.  4.  (a)  Intensity  distribution  produced  by  a  20-sphere-long 

chain  with  n  =  V3  as  a  result  of  illumination  with  collimated 
rays,  (b)  Linear  polarizer  (along  x')  placed  between  the  end- 
sphere  and  the  detectors,  (c)  Intensity  distribution  with  the 
polarizer  installed  also  showing  positions  of  two  detectors 
(dashed  rectangles)  used  for  DOP  calculations. 


it  follows  that  the  measure  of  the  intensity  of  azimuthally 
polarized  beams  can  be  obtained  by  using  a  90°  rotated 
local  detector  labeled  Icp  in  Figs.  4(b)  and  4(c).  Using 
relationships  and  ly  =  in  Eq.  (^,  we  calcu¬ 

lated  the  degree  of  radial  polarization  as  a  function  of 
n  for  chains  with  different  lengths.  We  ensured  good  con¬ 
vergence  of  the  results  with  reduction  of  the  width  of  the 
detectors.  Numerical  modeling  was  performed  by  using 
ZEMAX-EE. 

The  results  of  calculations  of  DOPs  for  three  types  of 
ray  sources  are  represented  in  Figs.  5(a)-5(c)  in  direct 
correspondence  with  the  structures  illustrated  in 
Figs.  3(a)-3(c).  It  is  seen  that  DOP  increases  with  the 
length  of  the  chain  reaching  ^-0.9  for  10-sphere  chains 
with  1.68  <n  <  1.80.  The  highest  DOP  can  be  obtained 
for  collimated  input  beams,  as  shown  in  Fig.  5(a).  In  the 
case  of  diverging  input  beams,  illustrated  in  Fig.  5(c), 
the  DOP  values  are  slightly  reduced  compared  to  that 
in  Fig.  5(a). 

The  DOP  calculations  for  the  case  of  spherical  source 
of  rays  are  illustrated  in  Fig.  5(b).  The  spherical  source 
produces  a  variety  of  rays  with  different  directionality. 
A  significant  fraction  of  these  rays  leaks  out  of  the  chain 
in  the  first  few  spheres  adjacent  to  S-sphere,  as  shown 
in  Fig.  3(b).  It  requires  longer  propagation  distance  in  this 
case  to  establish  a  typical  ring  pattern  in  the  transverse 
intensity  distribution  similar  to  that  illustrated  in 
Fig.  4(a).  Such  ring  patterns  indicate  that  the  PFMs 
are  formed  far  away  from  the  S-sphere,  typically  for  N  > 
10  [1].  In  addition,  the  spherical  source  is  coupled  to  both 
fundamental  PFMs,  one  of  which  is  shifted  by  the  sphere 
diameter  D  [1].  For  this  reason,  the  typical  transverse 
intensity  distribution  includes  the  ring  representing 
one  fundamental  PFM  and  the  central  intensity  peak 


Fig.  5.  (a)-(c)  Degree  of  radial  polarization  versus  n  calcu¬ 

lated  for  three  types  of  ray  sources  illustrated  in  Figs.  3(a)-3(c), 
respectively. 
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representing  the  second  fundamental  PFM  shifted  by  D. 
Formation  of  the  ring  also  makes  it  possible  to  apply  our 
technique  of  DOP  determination  based  on  using  the 
polarizer  with  two  stripe-detectors.  In  the  case  of  a 
spherical  source,  the  DOP  values  for  a  10-sphere  long 
chain  are  found  to  be  smaller  compared  to  that  for  other 
sources,  as  seen  in  Fig.  5(b).  However,  obtaining  DOP  in 
excess  of  0.9  is  still  possible  in  this  case  by  using  longer 
chains  with  N  =  20. 

In  conclusion,  we  showed  that  chains  of  cylinders 
and  spheres  have  a  common  property  of  low-loss,  high- 
extinction  ratio  polarization  filtering  of  PFMs  with  a  2D 
period.  Previously,  the  lossless  propagation  was  demon¬ 
strated  for  rays  incident  at  Brewster  angles  on  spheres 
with  n  =  Vs  ^  1.73  [1].  In  this  work,  we  showed  that  this 
effect  should  be  present  for  a  range  of  sphere  indices 
1.68  <?^<  1.80. We  showed  that  the  chains  of  cylinders 
filter  linearly  polarized  beams,  whereas  chains  of  spheres 
filter  radially  polarized  beams.  The  optimal  refractive 
index  for  observation  of  these  effects  is  very  close  to 
the  index  of  one  of  the  best  quality  infrared  materials, 
sapphire  (or  ruby)  with  n  ^  1.71  at  A^S  fim,  which  means 
that  various  practical  IR  focusing  and  polarization  devi¬ 
ces  can  be  built  using  sapphire  cylinders  or  spheres.  In 
addition,  since  focusing  of  radially  polarized  beams  can 
be  sharper  than  the  focusing  of  unpolarized  or  linearly 
polarized  beams  [1^,1^],  such  purely  dielectric  structures 
can  be  very  useful  in  high-resolution  polarization  imaging 
and  in  the  design  of  various  devices  including  focusing 
microprobes,  laser  scalpels,  and  optical  tweezers. 
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